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Although the 77-meson was discovered 40 years ago, only recently particle and nuclear 
physicists focused their attention on it. In many respects the //-meson is similar to the 
7T°-meson despite it being four times heavier. Both are neutral, spinless, and have almost 
the same lifetime, ~ 10 -18 sec. The kinship between the two mesons manifests itself very 
clearly in their decay modes. They are the only mesons which have a high probability of 
pure radiative decay, i.e., their quarks can annihilate into on-shell photons. The pion almost 
entirely decays into the radiative channel 7r° — > 7 + 7 (98.798%). For the 77 the purely 
radiative decay is also the most probable mode fl[], 
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7 + 7 (38.8%) 

7r° + 7r° + 7r (31.9%) 

7r+ + n- + n° (23.6%) 

7T+ + 7T- + 7 ( 4.9%) 

other decays ( 0.8%) 



> 

Therefore, when ir and r\ are viewed as elementary particles, they look quite similar. How- 
ever when one considers their interaction with nucleons, their difference is clearly manifested. 
Firstly, one expects a manifestation of the large ?77r -mass difference in the meson-nucleon 
dynamics and, at low energies, this is indeed observed. For example, the Sn-resonance 
iV*(1535) is formed in both ttN and 7]N systems, but at different collision energies, 

O 

E%fi{S n ) = 1535 MeV -m N -m n « 458 MeV 



E${S n ) = 1535 MeV -m N -m v « 49 MeV. 

Note that due to the large mass of the ?7-meson (547.45 MeV), this resonance is very close to 
the r/iV-threshold. Furthermore it is very broad, with V « 150 MeV, covering the whole low 
energy region of the 7]N interaction. As a result the interaction of nucleons with 77-mesons 
in this region, where the S'-wave interaction dominates, is much stronger than with pions. 
Another consequence of the Sn dominance is that the interaction of the 77-meson with a 
nucleon can be considered as a series of formations and decays of this resonance as shown 
in Fig. 1. 
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Fig. 1 

As with any resonant state, the A r *(1535)-resonance has branching ratios of the decay modes 
which do not depend on the formation channel and after its creation it decays into r)N and 
ttN channels with equally high probabilities JIJ 

( N + rj (35 - 55 %) 

A^*(1535) ->• | N + n (35 - 55 %) (1) 

[ other decays (< 10 %) 

Therefore, the series depicted in Fig. 1, must also include terms describing real and virtual 
transitions into the 7rA^-channel (see Fig 2). 




Fig. 2 

Thus in the energy region covered by the Sii -resonance, the r)N and irN interactions should 
be treated as a coupled channel problem. When such an analysis was performed, it was 
found that the near-threshold r)N interaction is attractive 0. This raises the question as 
to whether this attraction is strong enough so that an //-mesic nucleus can be formed. 

Since //-mesons decay very fast, it is impossible to produce beams of them and therefore 
they can only be observed in final states of certain nuclear reactions with other particles. 
This makes investigations of ?7-meson dynamics quite complicated. Therefore the possibility 
of sustaining an 77-meson inside a nucleus would be an exciting one as it would expose itself 
for a relatively long period in a series of successive interactions with nucleons, i.e., inside 
the nucleus it would undergo a series of absorptions and emissions through formations and 
decays of the A r *(1535)-resonance as depicted in Fig. 3. 




Fig. 3 
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The lifetime of such an ^-mesic nucleus would not be limited by the lifetime of the 
77-meson itself because after each creation of the Sn-resonance the ?7-meson is generated 
anew. However, such an 77-nucleus state can not be stable, since eventually the iV*(1535)- 
resonance will produce a pion instead of r\ as their creation probabilities are, according to 
(0), equally high. Of course such a pion can generate an iV*(1535)-resonance again which 
in turn may revive the 77 but such a possibility is rather low since the pion acquires, through 
the decay of the resonance, a kinetic energy of ~ 400 MeV and can thus easily escape. It is 
therefore clear that if an ^-meson is bound inside a nucleus, it can only be in a quasi-bound 
state with nonzero width. 

First estimation, obtained in the framework of the optical potential theory, put a lower 
bound on the number of nucleons A which is necessary to bind the ?7-meson, namely, A > 12 
0. Thereafter other theoretical investigations were devoted to this problem. All of them 
predicted 7/-nucleus bound states obeying this constraint. However, the search for narrow 
7/-nuclear bound states in an experiment with lithium, carbon, oxygen, and aluminum by 
Chrien et al. [f| produced negative results. 

The conclusion of this experimental work, however, did not discourage theoreticians in 
examining the possibility of an //-nucleus formation. The relatively large scattering lengths 
obtained for r/ 3 He and ?? 4 He systems using a zero-range r/iV-interaction [BJ cast doubt on 
the A > 12 constraint. Speculations of this kind are based on the argument that in the 
vicinity of the origin of the complex momentum plane the amplitude / can be replaced by 
the scattering length a and therefore the S"-matrix in this area can be written as 

1 + ika 

l + 2ikf^ 
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1 — ika 

This expression is valid only for small k and can have a pole in this region only if a is large. 
If a is negative the pole would be on the positive imaginary axis (bound state). Thus, a 
large negative scattering length indicates that a weakly bound state exists. Decreasing the 
interaction strength transforms the bound state into a resonance, and vice versa, implying 
a change of the scattering length from —00 to +00. 

Such simple reasoning, however, is valid only when the interaction is described by a real 
potential. In the case of an ^-nucleus system the inelastic rjA — > nA channel is always 
open giving rise to a significant imaginary part in the ^-nucleus potential. The resonance 
and quasi-bound state poles of the 5-matrix generated by a complex potential have quite 
different distribution in the complex fc-plane. In Ref. || it was shown that starting from a 
purely real potential and introducing an imaginary part which is gradually increased, results 
in S'-matrix pole-behaviour shown in Fig. 4. 
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Fig. 4 



Fig. 5 
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Therefore, in the case of a complex potential, both resonance and quasi-bound state poles 
are situated in the second quadrant of the complex momentum plane, under and above 
its diagonal respectively. The diagonal separates them because the energy E = k^/2fi 
corresponding to a pole at k = k , 



E ° = Y~ [( Rek o) 2 - (lmk ) 2 + 2i(Rek ){lmk ) 



2fi 

has a positive (negative) real part when k is under (above) it. Thus, the transition from 
resonances to quasi-bound states is a crossing of the diagonal. Since this can take place 
rather far from the point k = 0, we should not expect, in contrast to the real potential case, 
to have a large scattering length even if the binding energy, |Rei?o|, is small. Moreover, 
crossing the diagonal is not associated with dramatic changes of a. In short, scattering 
length calculations cannot provide a definite answer and a more rigorous approach must 
be employed. The most adequate way to solve this problem is to locate the poles of the 
^-matrix in the second quadrant of the fc-plane. In Refs. we developed a microscopic 
method that enabled us to calculate the elastic scattering amplitude for any complex value 
of k and thereby to locate its poles. The influence of inelastic channels is taken into account 
via a complex r]N potential. In what follows this method is described in somewhat more 
detail. 

Consider the scattering of an //-meson from a nucleus consisting of A nucleons. The 
Hamiltonian of the system is given by 

H = H + V vA + H A (2) 

where Hq is the free Hamiltonian corresponding to the 7/-nucleus relative motion, V v a = 

V\ + V2 + h Va. is the sum of rjN potentials, V{ = V v n(\R — f£|), where R and r*j are the 

coordinates of the rj and the z-th nucleon with respect to the cm. of the nucleus, and Ha is 
the total Hamiltonian of the nucleus, 

H A = -it V* + £ V NN (\ri ~ rj\) ■ (3) 

The elastic scattering amplitude f(k', k; z) describing the transition from the initial, \k, i/jq), 
to the final, \k',ipo), asymptotic state where \ipo) is the nuclear ground state and k the 
//-nucleus relative momentum, can be expressed in terms of the T-matrix elements 

f(k',k;z) = <k',ip \T(z)\k,ip > . (4) 

The operator T is related to the Green function Ga[z) = (z — Hq — Ha)^ 1 by 

T(z)=V + VG A (z)T(z) . (5) 

The task of solving Eq. (|5]) is a formidable one and thus one must resort to approximations. 
One such approximation is the so-called Finite-Rank Approximation (FRA) of the Hamil- 
tonian. It has been proposed in Refs. PJlOU as an alternative to the multiple scattering and 



optical potential theories. In this method the auxiliary operator 
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T°(z) = V + VG (z)T°(z), (6) 

where Gq(z) = (z — Hq)' 1 is the free Green function, is introduced. Using the identity 
A' 1 - B- 1 = B- l (B - A) A' 1 with A = z - H - H A and B = z — H , one gets the resolvent 
equation 

G A (z) = G {z) + G (z)H A G A (z) , (7) 

and thus 

T(z) = T°(z) + T (z)G (z)H A G A (z)T(z) . (8) 
The latter equation has the advantage that the spectral decomposition for the Hamiltonian, 

iJ A = ^£ n |i><i|+/^%Xfe|, (9) 

can be employed to bring Eq. (^j) into a manageable form. The FRA method is based on 
the approximation 

H A « £ |V>o >< ^o|, (10) 

which means that during the scattering of the ?7-meson, the nucleus remains in its ground 
state \ipo >• Such an approximation is widely used in multiple scattering and optical poten- 
tial theories where is known as the coherent approximation. Using (10) we get 

TO) = T\z) + £ T°(zMo > G (z)G (z - £ ) < ip \T{z) . (11) 

The matrix elements T(k', k; z) =< k' ,i(; Q \T(z)\k,il ! o > are thus given by 

T(k',k; z) = <k',ij \T°(z)\k,^ > 

+ S ° J W¥ (z-k>>>/2,)(z-£ -k>>>/2,) T{k ^ Z) - (12) 

The auxiliary operator T° describes the scattering of the ^-meson from nucleons fixed in their 
space position within the nucleus. This is clear since Eq. @ does not contain any operator 
acting on the internal nuclear Jacobi coordinates denoted by {r} = {xi,x 2 , - ■ ■ ,x A _i}. 
Therefore all operators in Eq. (||) are diagonal in the configuration subspace {r} and thus 

Tfi, E; r; ,) = V(k', k; f) + / T°(*>; k; f; z) (13) 

where 

< k', f'\T°(z)\k, f>= 8(f' - f) T°(k', k; r;z), < k', f'\V\k, r >= 5(r' - r) V(k', k; r) . 

It is clear that T°(k', k; r; z) depends parametrically on {r}. Therefore the matrix elements 
< k', ipo\T°(z)\k, ipo > can be obtained by integrating over the Jacobi coordinates 
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< k'^ \T°(z)\k^ >= / df]M^\ 2 T°(k',k;r;z). (14) 



Thus the solution of the scattering problem can be obtained by solving first Eq. (|T3|). 
averaging as in Eq. (|i~4|), and finally calculating T from Eq. (|i~2|). We must emphasize that 
the above scheme, is not the same as the first order optical potential approach used in the 
traditional pion- nucleus multiple scattering theory fill . Indeed, the latter is based on three 
approximations; i) the Impulse Approximation; ii) the omission of higher order rescattering 
terms in constructing the optical potential, and iii) the coherent approximation. In contrast, 
in the scheme considered here, the Impulse Approximation to obtain the rjN amplitude in 
nuclear media is not needed and no rescattering terms are omitted. 

The parameter z in the above equations corresponds to the total ^-nucleus energy, z = 
E — \Sq\ + iO, where E is the energy associated with the ^-nucleus relative-motion. On the 
energy shell we have E = k 2 /2fi. Therefore, even the auxiliary T°-matrix differs from the 
conventional fixed-scatterer amplitude in that it is always taken off the energy shell. In the 
case of scattering length calculations, we have E = and thus z = — \Sq\. This makes Eqs. 
(|j~2|) and ([13]) nonsingular and easy to handle. 

For practical calculations we rewrite Eq. (p~3| ) using the Faddeev-type decomposition 

T°(k',k;r;z)=J2T^k',k;r;z), 
i=i 

T^k',k;f;z)=U(k',k;f;z) + / ( ^)3 ^_ Effff. k; f; z) , (15) 

where tj is the t-matrix for the 77-meson scatterred by the nucleon i and is expressed in terms 
of the two-body t^-matrix via 

ti(k' , k;r; z) = t V N(k' , k; z) exp i{k — k') ■ r*j . (16) 

Expanding the \k,f > basis in partial waves and using the fact that, at the low energies 
considered here, the i]N interactions is dominated by the Sn-resonance, we may retain the 
S-wave only. The total orbital momentum is zero and since the ?7-meson is a spinless particle 
the nuclear spin can be ignored. Therefore, when Eq. (|1|) is projected on the S-wave basis 
\k,r >, it reduces to 



Tf(k>, k; r; z) = U{k\ k; r;z) + — jf dk" - ^ ^/^ ' E k; r, ; z) (17) 

where 

< k',r'\T°(z)\k,r >= 5 -^P- T?(k',k;r;z) 

and similarly for < k',r'\ti(z)\k,r >. 

The above formulae are given for the general case of A nucleons. In what follows we 
restrict ourselves to A =2, 3, and 4. The relevant Jacobi vectors are shown in Fig. 5. 



According to Eq. (0), U depends on the space configuration of the nucleons because k and 



k' are the 77-meson momenta with respect to the nuclear centre of mass while the nucleon i 
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is shifted from it by the vector r*j = a^xi + biX~2 + c%x% , where a% = |, 02 = — |, &i = 62 — 
Ci = c 2 = for the deuteron case; a\ = |, a 2 = — |, a 3 = 0, 61 = 6 2 = §, ^3 = _ 1, c i = C2 = 
c 3 = for the three-nucleon case, and ax = ~, a 2 = — ~, 03 = 04 = 0, b\ = b 2 = |, 63 = 64 = 
— ~, Ci = c 2 = 0, c 3 = |, c 4 = — I for the four-nucleon case. 
The S-wave projection of Eq. ( [16]) gives 

r dk'dk- - - 

(k',r I U(z)\k,r) = J ^ * dr" dr'" (k' , r ")(k' h f"\U(z)\ki, r'") (fa, r'"\k, r) 

Sir — t) 

= 2 3o( a i k ' x i)3o( b ik'x2)jo(cik'x 3 ) t vN (k', fc; z) j (a i kx 1 )j (b i kx 2 )jo(cikx 3 ) , 

where jo is the spherical Bessel function. The r)N interaction can be described by the 
t-matrix 

k] Z) = (k' 2 + a 2 )(z-E + tT/2)(k 2 + a 2 ) ' (18) 

This ansatz is motivated by the S^-resonance dominance. The vertex function for <-> N* 
is chosen as l/{k 2 + a 2 ) which in configuration space has a Yukawa-type behaviour. The 
propagator is taken to be of a simple Breit-Wigner form. With such a choice the ti has the 
following separable form 

U(k', k, r; z) = Hi{k'] r) t(z) H^k, r) (19) 

where 

, v A iQ(a i kx 1 )iQ{b i kx2)jo{.c i kxz) 

t{z) = — — — - , H i {k,r) = . 20 

2 — E + iF/2 (k z + cr) 

Therefore, 

T°(k', k, r;z)=J2 H t (k'; r) A^z) H 3 (k, r) (21) 

where 

(A- 1 ),, = -g- - (1 - ^ )r, ; (r.:) , r<,(r,s) = ^ ^ dfc ^ Jft(fc;r)ffj(fc,r) . 

(22) 

Formally, the the last integral involves products of six Bessel functions. However, several 
of the coefficients eij, 6j, and q are always zero and therefore only products of at most four 
Bessel functions can appear in the expression for ry(r, z) with the required integrals having 
the general form 

7 (p,«, Vl «,) = y o dk (k2 + a2nk2 _ p2 _ ioy (23) 
To calculate the latter integral we introduce the auxiliary one 
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l(p,u,v,v,,6) = ] Q dk {k2 + Q2)2{k2 _^_. 0)ik2 + s2)w2 (24) 



and thereafter evaluate the limit 



7(p, it, v , w) — hmj(p, u, v, w, S) . (25) 

8— >0 

The result thus obtained is 

7(p, w, u, io ) = — - — + w + 2w) - 2g(u + v) + g(u + v - 2w) 
lbuvw z 

-g(u -v + 2w) + 2g(u -v)-g(u — v- 2w)} , (26) 

where 

9( s ) = ( p2 + a 2)2 \ si S n ( Im P) ^3 ex P [«P|s|sign(Im p)] 

_ »exp(-a| 8 |) 2 + + 2 + q2 _ Mp^+a^j 

2cr " crp z J 

with 

. , . f +1, for a > > v 

sign (a) = i ' r n (28) 

v ' I -1, for a < . v ' 

To obtain the necessary nuclear wave functions ?/>o we employed the Malfliet-Tjon NN- 
potential [12] and the integro-differential equation approach (IDEA) [[13,14] which, for S- 



wave projected potentials, is equivalent to the exact Faddeev equations. 

Using the above formalism, the position and movement of poles of the ^-meson-light nu- 
clei ( 2 H, 3 H, 3 He, and 4 He) elastic scattering amplitude in the complex fc-plane are studied. 
The two-body t-matrix is assumed to be of the form (|TJ) with E = 1535 MeV — {m N + m v ) 
and T = 150 MeV 0. The range parameter used, a = 2.357 fm -1 , was obtained in a 
two-channel fit to the ttN —>■ nN and vrA^ — ► rjN experimental data 0. The parameter A is 
chosen to provide the correct zero-energy on-shell limit, i.e., to reproduce the 7]N scattering 
length djjN, 

. 2ix 
^(0,0,0) = a nN . 

The scattering length a V N, however, is not accurately known. Different analyses [T5[] provided 
values for the real part in the range Rea^ € [0.27, 0.98] fm and for the imaginary part 
Ima^Ar G [0.19, 0.37] fm. Therefore in a search for bound states one must use values of a V N 
within these ranges. To achieve this we used = (CO-55 + z0.30) fm and vary £ until a 
bound state appears. 

The poles found with a n N = (0.55 + i0.30) fm are shown in Fig. 6. The corresponding 
energies and widths are given in Table 1. When Rea^ increases, all the poles move up and 
to the right, and when a resonance pole crosses the diagonal it becomes a quasi-bound pole. 
The minimal values of Re a n N which generate 'zero-binding' (the poles just on the diagonal) 
are given in Table 2. 

All these values are within the uncertainty interval Kea V N G [0.27,0.98] fm. Thus even 
the possibility of an rjd binding cannot be at present excluded. Most recent estimates of 
Rea^ I3 are concentrated around the value Rea,^ ~ 0.7 fm, which enhances our belief 



that at least the a-particle can entrap an ^-meson. 
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Re p (fm 



Imp(fm ) 




system 


(MeV) 


T (MeV) 


7? 2 H 


31.46 


59.38 


77 3 H 


10.91 


22.68 


77 3 He 


10.14 


22.70 


77 4 He 


-2.05 


7.48 



Fig. 6 



Table 1 



For each of the four nuclei considered, the scattering lengths were calculated with eight 
values of the strength parameter A corresponding to Rea^: {(0.2 + O.ln) fm; n = 1,8}, 
which extends over the uncertainty interval. The Ima^ was fixed to the value 0.3 fm. 
An increase of Rea^jv moves the points along the trajectories, shown Figs. 7 and 8, anti- 
clockwise. When Rea^Ar exceeds the critical values given in Table 2, the rjN interaction 
becomes strong enough to generate a quasi-bound state. The corresponding 7/-nucleus 
scattering lengths are shown by filled circles (the trajectories for 3 He and 3 H are practically 
the same). 
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-4-3-2-1 IRe a vA 
Fig. 7 
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system 


minlRea^Ar} (fm) 


7? 2 H 


0.91 


r/ 3 H 


0.75 


?7 3 He 


0.73 


7/ 4 He 


0.47 



-4-3-2-1 lRea^ 
Fig. 8 



Table 2 



Finally, we would like to emphasize that the spectral properties of Hermitian and non- 
Hermitian Hamiltonians are quite different. Locating quasi-bound states is a delicate prob- 
lem which can be treated only by rigorous methods. As we have shown in Ref. || the r/A 
scattering length tells us nothing about the existence or not of an 77A quasi-bound state. 
This is clearly seen on Figs. 7 and 8 where the trajectories go smoothly from open to filled 
circles without any drastic changes or extreme values. 

In summary, it is shown that within the existing uncertainties of the elementary 77N 
iteraction all light nuclei considered can support a quasi-bound state which can result in an 
7/-mesic nucleus which is analogous to hypernuclei. Due to the specific quantum numbers 
of the ?7-meson (1=0, S=0) such states, if they do exist, can be used to access new nuclear 
states inaccessible by other mesons such as pions and kaons. Furthermore it can be used to 
elucidate the role played by the rj meson in Charge Symmetry Breaking reactions and in the 
violation of the Okubo-Zweig-Iizuka (OZI) rule. 
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